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Abstract 

We present a formulation of AA = (1, 1) super Yang-Mills theory in 
1+1 dimensions at finite temperature. The partition function is con- 
structed by finding a numerical approximation to the entire spectrum. 
We solve numerically for the spectrum using Supersymmetric Discrete 
Light-Cone Quantization (SDLCQ) in the large- A'c approximation and 
calculate the density of states. We find that the density of states grows 
exponentially and the theory has a Hagedorn temperature, which we 
extract. We find that the Hagedorn temperature at infinite resolution 
is slightly less than one in units of ^Jg'^Nc/ir. We use the density 
of states to also calculate a standard set of thermodynamic functions 
below the Hagedorn temperature. In this temperature range, we find 
that the thermodynamics is dominated by the massless states of the 
theory. 



1 Introduction 



In A/" = 4 super Yang-Mills theory at large there is a known mismatch be- 
tween weak and strong coupling by a factor of 3/4 in the free energy [1] . The 
weak-coupling result is calculable in perturbation theory; the strong-coupling 
result comes from black-hole thermodynamics. It would be interesting to be 
able to directly solve this theory at all couplings and see the transition be- 
tween the weak-coupling and strong-coupling limits Analytically this is 
generally not possible, although there have been a number of early discussions 
of methods for finite-temperature solutions to supersymmetric quantum field 
theory [3] . We will instead consider a numerical method based on Supersym- 
metric Discrete Light-Cone Quantization (SDLCQ) [H [S], which preserves 
the supersymmetry exactly. Currently this is the only method available for 
numerically solving strongly coupled super Yang-Mills theories. Conven- 
tional lattice methods have difficulty with supersymmetric theories because 
of the asymmetric way that fermions and bosons are treated, and progress [B] 
in supersymmetric lattice gauge theory has been relatively slow. 

Given that SDLCQ makes use of light-cone coordinates, with = {x^ + 
x^)/\/2 the time variable and p~ = (p" — p^)/\/2 the energy, we must take 
some care in defining thermodynamic quantities. The seemingly natural 
choice [7] of e~^LcP as the partition function has been shown by Alves and 
Das [8] to lead to singular results for well known quantities that are finite 
in the equal-time approach. They argue that using e~^^cP for the partition 
function implies that the physical system is in contact with a heat bath that 
has been boosted to the light-cone frame and that this is not equivalent to 
the physics of a system in contact with a heat bath at rest. 

A more direct way to see this is that, since the light-cone momentum 
p+ = [p^ _)_p3^^y2 is conserved, the partition function must include it in 
the form Z = e~^Lc(p with fi its chemical potential. The interpreta- 

tion of the chemical potential is that of a rotation of the quantization axis. 
Thus yU = 1 corresponds to quantization in an equal-time frame, where the 
heat bath is at rest and the inverse temperature P = \/2P^^, and /i 7^ 1 
corresponds to quantization in a boosted frame where the heat bath is not 
at rest. Thus /i corresponds to a continuous rotation of the axis of quanti- 
zation, and /i = would correspond to rotation all the way to the light-cone 
frame. A rotation from an equal-time frame to the light-cone frame is not a 
Lorentz transformation. It is known that such a transformation can give rise 
to singular results for physical quantities. This appears to be consistent with 



1 



the results found in [S]. A number of related issues have been extensively 
discussed by Weldon [9]. The method has also recently been apphed to the 
Nambu-Jona-Lasinio model [TO] . 

These difficulties are avoided if we compute the equal-time partition func- 
tion Z = e~^^ , as was proposed much earlier by Elser and Kalloniatis pT] . 
The computation may, of course, still use light-cone coordinates. Elser and 
Kalloniatis did this with ordinary DLCQ [121 US] as the numerical approxi- 
mation to (l-l-l)-dimensional quantum electrodynamics. Here we will follow 
a similar approach using SDLCQ to calculate the spectrum of A/" = (1, 1) 
super Yang-Mills theory in 1 + 1 dimensions [14J. Though this calculation 
is done in 1+1 dimensions, it is known that SDLCQ can be extended in a 
straightforward manner to higher dimensions ^TEl \TE\ [T7] . 

We have discussed the SDLCQ numerical method in a number of other 
places, and we will not present a detailed discussion of the method here; for a 
review, see |5]. For those famihar with DLCQ [121 US], it suffices to say that 
SDLCQ is similar; both have discrete momenta and cutoffs in momentum 
space, x~ G [—L,L]. In 1+1 dimensions the discretization is specified by a 
single integer K = (L/7r)P^, the resolution, such that longitudinal momen- 
tum fractions are integer multiples of 1/K. However, SDLCQ is formulated 
in such a way that the theory is also exactly supersymmetric. Exact super- 
symmetry brings a number of very important numerical advantages to the 
method; in particular, theories with enough supersymmetry are finite. We 
have also seen greatly improved numerical convergence in this approach. 

In Sec. [2] we briefiy review super Yang-Mills theory in 1+1 dimensions. 
The discussion in Sec. |3]describes the method we use to extract the density of 
states from the numerical spectrum. The calculation of the density of states 
is presented in Sec. HI We fit the data to smooth analytical functions, and we 
find that the theory has a Hagedorn temperature Th [LSI, which we calculate. 
In Sec. ISl we use the analytic fit to the density of states to calculate the free 
energy, the energy, and the specific heat, up to the Hagedorn temperature. 
Section [6] contains a discussion of our results and the prospects for future 
work using these methods. 
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2 Super Yang— Mills theory 



We will start by providing a brief review of M = (1, 1) supersymmetric 
Yang-Mills in 1+1 dimensions. The Lagrangian of this theory is 



C = Tt 



(1) 



The two components of the spinor \Ef = 2~^/^(^) are in the adjoint represen- 
tation, and we will work in the large- iVc limit. The field strength and the co- 
variant derivative are F^i. = d^Ay-d„A^ + ig[Afj,, A„] and = dfj, + ig[A^, ]. 
The most straightforward way to formulate the theory in 1+1 dimensions is 
to start with the theory in 2+1 dimensions and then simply dimensionally 
reduce to 1+1 dimensions by setting (j) = A2 and 82 ^ for all fields. In the 
light-cone gauge, A'^ = 0, we find 



(2) 



The mode expansions in two dimensions are 



0ii(O,X ) 



27r Jo V2fc+ 



1 



dk^ 



hij{k' 



+ hUk 



2^ Jq 

To obtain the spectrum, we solve the mass eigenvalue problem 

2P+P-|V9) = V2P+{Q-f\^) = M^\^). 



(3) 



(4) 



This theory has two discrete symmetries, besides supersymmetry, that 
we use to reduce the size of the Hamiltonian matrix we have to calculate. 
S'-symmetry, which is associated with the orientation of the large-iV^ string 
of partons in a state [TH], gives a sign when the color indices are permuted 



(5) 



S : aij{k) -> -aji{k), bij{k) -hji{k). 



P-symmetry is what remains of parity in the direction after dimensional 
reduction 

P : aijik) -aij{k), hjik) hjik). (6) 
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All of our states can be labeled by the P and 5* sector in which they appear. 

We convert the mass eie: envalue problem 2P+p-\M) = M'^\M) to a ma- 
trix eigenvalue problem by introducing a discrete basis where P~^ is diagonal. 
We will always state in units of g'^Nc/ir. In SDLCQ the discrete basis is 
introduced by first discretizing the supercharge Q~ and then constructing P~ 
from the square of the supercharge: P^ = {Q^Y / \/2. To discretize the the- 
ory, we impose periodic boundary conditions on the boson and fermion fields 
alike and obtain an expansion of the fields with discrete momentum modes. 
We define the discrete longitudinal momenta as fractions nP^/K of the 
total longitudinal momentum P"*", where K is an integer that determines the 
resolution of the discretization and is known in DLCQ as the harmonic reso- 
lution [12]. Because light-cone longitudinal momenta are always positive, K 
and each n are positive integers; the number of constituents is then bounded 
by K. The continuum limit is recovered by taking the limit K —>■ oo. 

In constructing the discrete approximation, we drop the longitudinal zero- 
momentum mode. For some discussion of dynamical and constrained zero 
modes, see the review [I3] and previous work [20]. Inclusion of these modes 
would be ideal, but the techniques required to include them in a numeri- 
cal calculation have proved to be difficult to develop, particularly because 
of nonlinearities. For DLCQ calculations that can be compared with ex- 
act solutions, the exclusion of zero modes does not affect the massive spec- 
trum [13]. In scalar theories it has been known for some time that con- 
strained zero modes can give rise to dynamical symmetry breaking [T^ , and 
work continues on the role of zero modes and near zero modes in these the- 
ories [211 [221 [231 [22 [25]. 

3 Density of states 

The thermodynamic functions will be written as sums over the spectrum 
{Mn} of the theory. The most convenient way to calculate such sums is to 
represent each sum as an integral over a density of states p(M^), 



From our numerical solutions we can approximate the density of states by 
a continuous function. The remaining integrals in the thermodynamic func- 



n=l 




(7) 
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tions are then done by standard numerical integration techniques, which are 
fast and convenient. 

We can look at the density for a series of increasing resolutions K in the 
SDLCQ approximation and thereby discuss the convergence of the density 
in the limit K —>■ oo. The maximum mass that we can reach in the SDLCQ 
approximation increases as we increase the resolution. We report results for 
11<K< 17. 

Convenient functions to extract from the spectral data [26] are the cu- 
mulative distribution function (CDF) N{M'^,K) and the normalized cumu- 
lative distribution function (NCDF) /(M^, K, M^). The CDF is the number 
of massive states at or below at resolution K, and the NCDF is this 
number divided by the total number of massive states below an arbitrary 
normalization point M^, again at resolution K: 

The function / turns out to be very smooth and can be fit by a single smooth 
analytic form. By definition, the density of states is given by 

It is also convenient to define a normalized density of states [26] 

PiM^K,M^;)^'-l^Ml^. (10) 

It is well known that if the density of states grows exponentially with the 
mass of the state, 

p(M2)~exp(M/TH), (11) 

the theory will have a Hagedorn temperature, Th [18j. Above the tempera- 
ture Th, the thermodynamic integrals diverge. 



4 Numerical results for the density of states 

The numerical results presented in this section are the first from our new 
code, which was rewritten to run on clusters. Most of these results were 
produced on our six-processor development cluster. While this cluster was 



5 



sufficient for this problem, we expect to be able to handle larger problems by 
moving to larger clusters. In fact, it is now so easy to generate the Hamilto- 
nian up to resolution K = 17, that we only used one node in our cluster for 
that purpose. What made this calculation challenging numerically was that 
we needed to extract a large number of eigenvalues. For the largest values of 
the resolution, this was done with a specially tuned Lanczos diagonalization 
code based on the techniques of Cullum and Willoughby [2Z] • We should note 
that, prior to the development of the new matrix-generation code, the highest 
resolution results presented for this theory were for resolution K = 10 [13]. 
Here we will present only resolutions K > 11. All our earlier results are now 
trivially reproduced. 




50 LOO 150 200 250 300 



f (M^, K, 100) 
0.025 [ 




(b) 



Figure 1: The distribution functions (a) N{M'^,K) for K = 12 and (b) 
/(M2, K, 100) for K = 12, 14, and 17, as functions of in units of g^NjTr. 
In (a) we also show the error-function fit aerf(6(x — c)) + d, and in (b), a 
polynomial fit. 

In the SDLCQ approximation, the portion of the spectrum that we can 
see at any finite resolution will naturally be cut offi As we approach the 
cutoff, the approximation limits the number of states that are available and 
distorts the density of states. In Fig. [1^, we present the data for the CDF 
at resolution K = 12, where we can diagonalize the entire Hamiltonian. 
We clearly see evidence of this distortion. By the midpoint, the cutoff is 
already diminishing the number of states available in the approximation. 
Interestingly, we can find a fit to this data with a simple universal function 
of the form aerf(&(a; — c)) + d. Clearly the fit shown is excellent; the fit is 
so good that one cannot separately see the data and the fitted curve on this 
scale. 



6 



p(K, M^) 

0.04 r 



0.03 



0.02 



0.01 




20 



40 



60 



80 



100 




Figure 2: The normalized density of states (a) for < 100 and (b) extrap- 
olated to all masses. The data points are a convenient way of displaying the 
continuous functions calculated from the fits to the CDF. 



At low masses, there is a mass gap, which has been discussed elsewhere [5]. 
The mass gap closes linearly with 1/K. For very small values of M^, the 
decrease of the mass gap with increasing resolution adversely affects the 
quality of the universal fit, and it is convenient in that region to improve the 
quality of the fit by using a polynomial function of the form 

p{K) 

f{x,K, WO) = J2 apx''<+Pe{x-x^UK)). (12) 

p=0 

In Fig. [T)d we show the fit to the NCDF for some representative values of the 
resolution K. We have only shown the data at resolutions i^' = 12, 14, and 
17 to keep the figure uncluttered. One can see how the endpoints tend to 
lower values as the mass gap closes with increasing K. 

At larger values of K it is difficult to completely diagonalize the entire 
Hamiltonian. We have limited ourselves to states with < 100. However, 
once we know the universal form of the function that fits the CDF, we can 
fit just the region < 100 and extrapolate to all masses. At large M^, 
the CDF approaches the total number of bound states. The total number of 
states in the SDLCQ approximation at any resolution, and, in any symmetry 
sector in the large- Ai'c approximation, is exactly calculable; the general results 
will be discussed elsewhere. We use this asymptotic value of the CDF, in 
addition to the behavior for < 100, in making the fits. We have done 
this at all resolutions up to = 17. In Fig. [2^ we show the normalized 
density of states calculated from the NCDF, for < 100. In Fig. we 
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show the normahzed density of states for 11 < < 17, extrapolated to the 
full range of masses. 




Figure 3: The Hagedorn temperature Tu, as obtained from exponential fits, 
aexp — «, to the universal fits to N{M'^, K), as shown in (a) for K = 17, 
and (b) a linear plot against 1/K for fits in the range 11 < K < 17. 

Inspecting these curves, we see that on the up-slope part of the density 
of states, where we believe our numerical approximation is a valid approx- 
imation to the actual density of states, the shape appears exponential. As 
we go to larger and larger values of the resolution K, the size of this region 
grows. This suggests that the density of states ultimately becomes simply 
an exponential, and, therefore, this theory has a Hagedorn temperature. To 
find the Hagedorn temperature, we fit the NCDF in this up-slope region with 
a function of the form 

/(M2,if,100) = aexp(^^) (13) 

In Fig. [3b we plot Th against 1/K. This yields a good linear fit, which we 
extrapolate to infinite resolution. We find that the Hagedorn temperature at 
infinite resolution is slightly less than one in units of \J'(pN\Jn. This value 
serves as a limiting temperature for the region of validity in the calculation 
of thermodynamic quantities. 

5 Finite temperature in 1+1 dimensions 

In the large-iVc approximation, the numerical solution of a theory is a set 
of non-interacting bound states. Therefore, the thermodynamics of such 
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supersymmetric theories is simply the thermodynamics of a gas of a large 
number of species of degenerate bosons and fermions. In principle, one could 
go beyond the calculation of the standard set of the thermodynamic functions 
and calculate a variety of matrix elements. These calculations would require 
the wave functions of the bound states, which can be calculated as part of the 
SDLCQ calculation. We will, however, not exploit this detailed information 
here. We will focus on the calculation of standard thermodynamic quantities 
that can be obtained from the density of states. The light cone plays no 
role beyond the calculation of the density; the thermodynamics is that of a 
system at rest. 

Let us now briefly review the thermodynamics of free bosons and fermions. 

We assume that our system has constant volume V and is in contact with a 
heat bath of constant temperature T. The free energy in units with ks — 1 
is 

r{T,V) = -Tl-aZ. (14) 
The contribution of a single bosonic oscillator to the free energy Fg is 

Fb = 2VT r ^In f 1 - e-P°/^) , (15) 

where Pq = \J Af ^ + pi and the factor of 2 compensates for integrating over 
only positive values oip^. It is convenient to change variables from p^ to Pq: 

dpz = , ^° dpQ. (16) 

The limits of integration are changed from < < oo to M < po < oo- We 
may also use the following representation for the logarithm that appears in 
the integrand 

ln(l-e--o/^) = -E^^, (17) 

3=1 ^ 

since po is positive and < e^^"/^ < 1. Finally, we obtain an expression for 
the total bosonic free energy just by summing over the energy spectrum 

r,^.Ylttr,,„v^(^^\ (18) 

The calculation of the fermionic contribution to the free energy proceeds 
analogously, and we find the identical result with the exception that there is 
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a factor of (—I)-' inside the summation. We can separate out the massless 
states from these expressions and calculate their contribution explicitly. We 
know that for resolution K there are {K — 1) massless bosons and {K — 1) 
massless fermions. Thus the contribution to the free energy from massless 
states is 



-r-o 

B 



{K - l)7r 



{K - l)7r 



6 ' ^ 12 

After doing the integral over pq, we find for the total free energy 



V) {K ■ 



yrp2 



_l)7r_ Xi((2Z + l)f^) 

^ n=l 1=0 



{21 + 1) 



(19) 



(20) 



The even terms, where j = 21 in the original sum, cancel between the fermion 
and boson contributions. We have also factored out the temperature depen- 
dence of the massless contribution and the volume dependence. 
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Figure 4: The free energy (a) and the heat capacity (b) as functions of tem- 
perature for each resolution. Both functions are normalized to the number 
of massless states, Nq = 2{K — 1) . The data points are a convenient way to 
display the continuous functions calculated from fits to the CDF. 

We can now rewrite the free energy in terms of the density of states. 
The sums involving the Bessel function are cut off at a few terms; generally 
^cut < 2 will be sufficient. We find 



J^{T, V) _ {K-l)n _ J_ 



roc 



/cut Kl 

dM^p{M^)Mj2 — 

1=0 



{21 



l)f 



{21 + 1) 



(21) 
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The free energy may now be used to calculate all the thermodynamic 
functions. The internal energy and heat capacity are given by, 

It is straightforward, given the density of states, to calculate the thermo- 
dynamic functions. Fig. shows the free energy, and Fig. HJo shows the heat 
capacity. We expect the free energy to diverge as A*"^ and therefore must 
normalize our results to extract a finite number. In most of the region below 
the Hagedorn temperature the thermodynamic functions are totally domi- 
nated by the massless states. It therefore seems appropriate to normalize 
the thermodynamic functions to the total number of massless bound states, 
which is a function of the resolution and is 2{K — 1). Alternatively, we could 
normalize by the number of states in any region. It is conceivable that at very 
high resolutions, where the mass gap is significantly less than one, that the 
massive states may make an important contribution to the thermodynamics. 
In that case we would not choose to normalize by the massless states. 

6 Discussion 

The \axge-N(. SDLCQ solution of M = (1, 1) super Yang-Mills theory in 
1+1 dimensions gives a set of non- interacting bound states. From this set of 
bound states it is in principle possible to calculate the thermodynamics of this 
theory. Central to this calculation is the calculation of the density of states. 
At resolutions K = 12 and below, where we can completely diagonalize the 
Hamiltonian, we find that the entire cumulative distribution function can be 
fit with a single erf function. From the cumulative distribution function, it 
is straightforward to calculate the density of states. For K larger than 12, 
it is difficult to calculate the entire spectrum; therefore, our calculations are 
confined to a fixed range of masses, < 100 g'^Nc/ir. Using the known 
form of the distribution, we only need to fit a section of the cumulative 
distribution function to get a very good fit to the entire distribution. We 
know analytically the total number of bound states at any resolution, and 
this information can also be used in conjunction with a fit to a section of the 
distribution to produce the fit to the entire distribution. 

The density of states that are found by this procedure grow sharply at 
small masses, then level off and decrease at larger masses. The peak of the 
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density of states grows as we increase the resolution. Our understanding 
of this behavior is that the cutoff of the theory is forcing the density of 
states to level off, turn over, and then decrease. The true behavior of the 
density of states is reflected in the region of the density of states that is 
rapidly increasing, because it is this region that is increasing in size with the 
resolution K. 

It appears that the cumulative distribution function, and, therefore, the 
density of states, are growing exponentially with the mass. To confirm this 
and find the asymptotic values of this growth, we fit the cumulative distri- 
bution function with an exponential at each K. We extrapolate these results 
to infinite resolution to find the asymptotic behavior of the density of states. 
The coefficient of the exponential growth is the reciprocal of the Hagedorn 
temperature. We find that this temperature is about 0.854 \J g'^Nc/n/kB- 
The thermodynamic functions calculated from this data are expected to pro- 
duce valid results up to a temperature that is around Th- 

It is now straightforward to calculate a standard set of thermodynamic 
functions from this density of states. The best estimate of the thermodynam- 
ics is obtained by using the exponential fits to the density of states. What 
we see is that, for resolutions up to = 17, all of the massive bound states 
are well above the Hagedorn temperature. The thermodynamics below Th 
is therefore controlled by the K — 1 massless boson bound states and K — 1 
massless fermion bound states. 

We can speculate on what will happen as the resolution goes to infinity. 
We have seen that the mass gap closes linearly with 1/K. So, for a reso- 
lution of order 100, there will be massive bound states below the Hagedorn 
temperature. This, of course, assumes that the estimate of the Hagedorn 
temperature is not changed by the higher resolution calculations. We found, 
however, that the actual number of massive bound states in a fixed mass 
range may grow slowly. For resolutions 11 to 17 we are able to find excellent 
fits with both exponential and linear growth as a function of the resolution 
K for masses with < 100. If the number of massive states grows only 
linearly with K, the contribution to the thermodynamic functions below Th 
might become significant but not dominant. 

These calculations indicate that M = (1, 1) super Yang-Mills theory 
in 1 + 1 dimensions has a Hagedorn temperature of about one in units of 
\J g'^Nf./'n/kB- More generally, we found that SDLCQ can be used to find in- 
teresting properties of finite-temperature supersymmetric field theories. The 
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extension of this method to theories with more supersymmetry and in higher 
dimensions appears to be straightforward but may be computationally chal- 
lenging. 
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